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We measure the renormalized coupling in the Twisted Polyakov loop scheme for SU(3) gauge 
theory coupled with Nf = 12 fundamental fermions. To find the infrared fixed point of this theory, 
we focus on the step scaling function for the renormalized coupling. We take the continuum limit 
using the linear function of (a/L) 2 and a constant fit function. We find that there is a sizeable 
systematic error due to the choice of the continuum extrapolation function in the low energy 
region. We will give some directions to reduce the systematic errors. 
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1. Introduction 

Discovering the fixed point and the critical phenomena around it are one of the intriguing 
subjects in quantum field theory. The fixed point is defined by the zero point of the beta function of 
the theory, and there the theory exhibits scale invariance. Some of them are exactly solvable. The 
values of the coupling constant at the fixed point depend on the renormalization scheme. In general, 
changing the renormalization scheme corresponds to the coordinate transformation in theory space. 
Thus, the renormalization group flows are changed by this transformation, but the existence of fixed 
point and whose critical exponents are not changed[|l], g]. 

Recently, there have been many papers concerning the fixed point search and the study of 
phase structure of SU(N) gauge theory with a large number of flavors. The perturbative j8 function 
indicates the existence of non-trivial infrared (IR) fixed points for a certain range of the number of 
large-flavor (Nf) SU(N) gauge theories, which is so-called "conformal window". Possible appear- 
ance of these IR fixed points has stimulated phenomenological studies of topics such as dynamical 
electro-weak symmetry breaking and unparticle physics. 

The existence of these IR fixed points depends on the gauge group, the number of flavors, 
and the representation of fermion fields. For SU(3) gauge theory with fermions in the fundamental 
representation, such a fixed point has been predicted in the range 8 < Nf < 16 using perturbation 
theory [Q]. However, the value of the renormalized coupling may be in the regime where per- 
turbation theory is not applicable. There are some analytical studies with improvements of the 
perturbation for the search for the lower bound of the conformal window, and the largest one give 
a prediction that Nf < 12 is out of the window [Q]. 

One of the nonperturbative approaches to this subject is a lattice simulation, and there are 
some recent studies concerning it. First such lattice study for SU(3) gauge theory was carried out 
in Ref. [j5|, where the authors investigated the phase structure of the case of Nf = 16. Appelquist 
et al. performed lattice calculation of the running coupling constant in the Schrodinger functional 
(SF) scheme and discovered evidence of an IR fixed point in the case of Nf = 12 [§]. On the other 
hand, Fodor et al. does not obtain a signal of IR fixed point using the different scheme based on 
the Wilson loop[|7], §]. The difficulty is mainly due to scheme dependence of the running coupling 
constant and the presence of significant lattice artifacts in the strong-coupling regime. Therefore 
it is important to measure the running coupling in different renormalization schemes, and estimate 
the dicretization error carefully. 

In this work, we perform lattice simulation of the running coupling constant for the funda- 
mental representation, Nf = 12, SU(3) gauge theory. Similar to the approach of Appelquist et al., 
we measure the step scaling function a(s,g 2 (L)) = g 2 (sL) keeping the values of bare coupling 
constant (j3) that give constant renormalzed coupling (g 2 (L)) for each small lattice size. We work 
in the Twisted Polyakov Loop (TPL) scheme which does not contain 0(a/L) discretization errors. 
This TPL scheme was first proposed by de Divitiis et al. [|9|, [K)|] for SU(2) gauge theory, and we 
extend the definition of the scheme to the SU(3) case[fnj]. 

In this paper, we give a short review of TPL scheme in §.||. Our preliminary results for Nf = 12 
SU(3) gauge theory is reported in §.[| In the last section, we discuss some ways of improvements 
to reduce the discretization error and show the future direction to search for the IR fixed point. 
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2. Twisted Polyakov loop (TPL) scheme 

In this section, we present the definition of the Twisted Polykov Loop scheme in SU(3) gauge 
theory. This is an extension of the SU(2) case as discussed in Ref. [B]. To define the TPL scheme, 
we introduce twisted boundary condition for the link variables in x and y directions on the lattice: 

U ti {x + vL/a) = a v U^(x)al. (v=x,y) (2.1) 

Here, D. v are the twist matrices.The gauge transformation for link variables Up (r) — > A(r)Up (r)A^ (r + 



p.) and eq.(2.1) imply the gauge transformation at boundary with a twisted gauge matrix: 

A(r + PL/a) = Q v A(r)Qj;. (2.2) 

Because of this twisted boundary condition, the definition of Polyakov loops in the twisted direc- 
tions are modified, 

P x (y,Z,t) = Tr ([]JU x (x = j,y,z,t)]a x e a,ly l^ , (2.3) 

in order to satisfy gauge invariance and translation invariance, and similary for y-direction. 

The renormalized coupling in TPL scheme is defined by taking the ratio of Polykov loop 
correlators in the twisted (P r ) and the untwisted (P z ) directions: 

2 1 (L yiZ P JC (y,z,L/2 a )P < (0,0,0)t) 

8TP k(^ y P z (x,y,L/2a)P z (0A0Y)' 

At tree level, this ratio of Polyakov loops is proportional to the bare coupling. The proportionality 
factor k is obtained by analytically calculating the one-gluon-exchange diagram. To perform this 



analytic calculation, we choose the explicit form of the twist matrices [|12|], and in the case of SU(3) 
gauge group we found k ~ 0.03184 jTT|]. 

To introduce the fermions which satisfy both the twisted boundary condition and translation 
invariance on the lattice, we have to introduce additional "smell" degree of freedom and indentify 
the fermion field as a N c x N s matrix (y/^(x)), where N c and N s are the numbers of color and smell 
degrees of freedom respectively [|13|]. Then we impose the twisted boundary condition for fermion 
fields to be 

r a (x + vL/a) =e i "/ 3 Q a v b y£(n v )J a (2.5) 

for v = x,y directions. Here, the smell index can be considered as a "flavor" index, then the number 
of flavors should be a multiple of N s (= N c = 3). We use staggered fermion in our simulation. This 
contains four tastes for each flavour. This enables us to perform simulations with a multiple of 
12-flovors in this SU(3) gauge theory with twisted boundary condition. 

At the end of this section, we would like to remark on the center symmetry of this theory. The 
generator of the symmetry is given by z = exp(2nik/3), k = 0, 1,2 for SU (3) gauge theory. Al- 
though the Wilson gauge action is invariant under the following transformation for the link variable 
for each direction, 

Up(t,x)^z.Up(t,x), (2.6) 
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the fermion action is not invariant. Therefore the vacuum expectation values of operator at z = 1 
and at z = exp{±2ni/3) are different. In the simulation, we generated the gauge configurations 
at nontrivial vacua where the vacuum expectation value of Polyakov loop correlator in untwisted 
direction has a nontrivial phase. We also investigated the one-loop effective potential for each 
vacuum analytically, and found the nontrivial one is the true vacuum] 14]. 



3. Simulation detail 
3.1 Simulation set up 

The gauge configurations are generated by the Hybrid Monte Carlo algorithm, and we use the 
Wilson gauge and the staggered fermion action. To reduce large statistical fluctuation of the TPL 
coupling, as reported in Ref. [JTJ], we measure Polyakov loops at every Langivin step and perform a 
jackknife analysis with large bin size, typically of 0(1O 3 ). This enables us to evaluate the statistical 
error correctly. The simulations are carried out with several lattice sizes (L/a = 4, 6, 8, 10, 12, 16) at 
more than twenty j8 values in the range 4.5 < j8 < 20. We generate 50,000-100,000 configurations 
for each (j3,L/a) combination. 
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Figure 1: TPL renormalized coupling in the each j3 and L/a in Nt = 12. The right figure shows the detail 
behavior in low /3 region. 

Figure [I] shows the j3 dependence of the renormalized coupling in TPL scheme at various 
lattice sizes. The results can be fitted at each fixed lattice size to the interpolating function 

2 /ox C l . C 2 C 3 C 4 

8tp(P) = J + ^ + ^I + p M 

where C, are the fit parameters. 

The right figure in Fig. [I] shows the detailed behavor of low-beta region. We find that gj P 
increases monotonously with the increase of the lattice sizes within the whole range of /3 examined. 
This holds even at /3 = 4.5 This is in the contrast to the case of Schorodinger functional scheme 
where the renormalized coupling exhibits some "crossing" behavior 

3.2 The continuum extrapolation and the scaling function 

We investigate the growth rate of TPL coupling in the continuum limit within the step scaling 
method. The procedure of the step scaling is first to find a set of bare coupling constant (j8) for 
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each small lattice size, which gives an input value of renormalized coupling (g 2 (/3 , a /L) = u). Next, 
we measure the step scaling function L(u,s,a/sL) = g 2 (p,sL/a)\ g 2^ a /L)=u- Finally, we take the 
continuum limit and obtain the step scaling function a(s,u) in the continuum. 

o(s,u) = HmL(u,s,a/sL)\ i (L)=u . (3.2) 

In our study, we use L/a = 4,5,6,8 as a small lattice size of the step scaling, and estimate the 
coupling constant for L/a = 5 from interpolations at the fixed j8 using the interpolation fit results of 
the lattice sizes L/a = 4, 6, 8. The step scaling parameter is s = 2, thus a large lattice size of the step 
scaling is L/a = 8, 10, 12, 16. From now, we denote the step scaling function a(s = 2,u) = c(u). 

The total error of the step scaling function can be estimated by the sum of the statistical error 
from each data and the systematic errors. The systematic errors have two origins. One of them is 
included in each data value L(u,s,a/sL), in which the value of /3 that is tuned to giving an input u 
have the error. The other one comes from the continuum extrapolation. If we measure the running 
behavior of coupling constant, the systematic error which comes from j8 -tuning is accumulated. 
However, we can carry out each step-scaling procedure independently for a given u. We focus on 
the growth rate of the renormalized coupling for several values of u. And, if the running coupling 
constant reaches the fixed point, the growth rate o{u)/u should be 1. 

Now, we have to consider the statistical error and the systematic error which come from the 
continuum extrapolation. To estimate the systematic error, we take the continuum limit using both 
a constant fit and a linear function in (a/L) 2 . In the case of quenched QCD, we found that the 
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Figure 2: The continuum limit of g TP with s = 2. In the left figure, each line denotes a linear fit function of 
(a/L) 2 . The statistical error bars are of the same size of the symbols (less than 3%). The each input value 
of TPL coupling (u = g\ p (L/a)) is u = 2.36,2.04, 1.73, 1.41, 1.10,0.79,0.47 from top to bottom. The right 
figure shows the detailed behavior of the lowest energy scale in the left figure. The pink line denotes the value 
of the input renormalized coupling u = 2.36. The green and blue lines denote the linear extrapolation for four 
points (L/a = 8,10, 12, 16) and the constant exUapolation for three points (L/a = 10, 12, 16) respectively. 

coupling constant of the TPL scheme exhibits scaling behavior even at the small lattice sizes, as 



shown in Ref. [11]. In the case of Nf = 12, we show the continuum extrapolation and scaling 
behavior for each step scaling in Fig. ^ We found that the scaling behavior in the low energy 
region becomes worse. The right figure in Fig. show the detailed behavior of the lowest energy, 
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in which the input renormalized coupling is u = 2.36. We found that the data of L/a = 8 is far from 
the data of the other lattice size and there is a large scaling violation in the step from L/a = 4 to 
sL/a = 8. 



!7 ((()/'[( 





Figure 3: The growth rate of the scaling function using a linear extrapolation (the left figure) and a constant 
extrapolation (the right figure). 



The growth rate of the scaling function (a(u)/u) in the continuum limit is shown in Fig. ||[ in 
which we use a linear extrapolation of (a/L) 2 and a constant extrapolation. 

The left figure in Fig. || may show that the growth rate is decreasing for u > 2.2. However, 
the right one shows that there is no such signal. Actually, the difference between these two extrap- 
olations gives a systematic error of the renormalized coupling. At the low energy region, there is 
more than 6% systematic error. The large systematic error comes from the large scaling violation 
of a small lattice size, as we show in Fig. |2[ To reduce the systematic errors and to conclude the 
existence of IR fixed point, we have to do some improvements or carry out the simulation for larger 
lattice size. 



4. Discussion and future direction 

In this paper, we investigate the growth rate of TPL coupling in low energy region j8 ~ 5.0. 
The TPL coupling does not show the inversion of the order of lattice size even in the low j8 region. 
The lattice renormalized coupling has a discretization error and it depends on the renormalization 
scheme. To remove the discretization effects and to estimate the systematic error, we took the 
continuum limit using two ways: a linear extrapolation of (a/L) 2 and a constant extrapolation. The 
statistical error is less than 3% even in a low energy region, but the systematic error is more than 
6%. We cannot give a conclusive statement for the existence of the fixed point due to the large 
systematic error. 

To solve this difficulty and to obtain the conclusion of the fixed point search in Nf = 12, there 
would be several directions to take. One direction is to carry out the simulations on larger lattices. 



In Ref. [16] we carry out the step scaling using a set of large lattice size, L/a = 6,8, 10, 12 and 
sL/a = 9, 12, 15, 18, for SU(2) eight flavor case. In that study, the continuum extrapolation behaves 
nicely even in the low energy region. Second direction may be given by some improvements of the 
action to reduce the discretized effect in the low energy. Third one is to estimate the discretization 
error using the analytical calculation by lattice perturbation theory, and then to subtract this value 
of discretization error from the lattice raw data. That may give a better scaling behavior. 
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Finally, we would like to show our future direction. We have to measure the anomalous di- 
mension of the field around the fixed point, if there is a fixed point in the infrared region. The 
anomalous dimension of the operator on the fixed point is related with the conformal algebra of 
the conformal field theory, so that it should be scheme independent. Comparision of its values 
which are measured in several schemes gives a conclusive statement of the existence of the fixed 
point. Furthermore, from the phenomenological motivation we expect that the fermion composite 
operator has a large anomalous dimension and gives a origin of Higgs sector of standard model. 



We will report a new method of measurement of the anomalous dimension in Ref. [14] 
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